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Abstract. We prove the Lukacs characterization of the Wishart distribution o n non octonion symmetric 
cones of rank greater than 2. We lessen the smoothness assumptions in the iBobecka and W csolowski 
[2002i | version of the Lukacs theorem. The main tool is a new solution of the Olkin-Baker functional 
equation on symmetric cones, which was solved under assumption of continuity of respective functions. 
It was possible thanks to the use of Gleason's theorem. 



1. Introduction 



The iLukacs 1955| theorem is one of the most celebrated results in the area of characterizations of 
probability distributions. It states that if X and Y are independent, positive, non- degenerate random 
variables such that their sum and quotient are also independent then X and Y have gamma distributions 
with the same scale parameter. This theorem has man y g eneralizations. The mos t important in the 
multivariate setting were given in lOlkin and Rubinl 19621 ] and lCasalis and Letad 1996J, where the authors 
extended characterization to matrix and symmetric cones variate distributions, respectively. There is no 
unique way of defining the quotient of elements of the cone and in these papers the authors have considered 
very general form U — [w(X + Y)] -1 X[w T (X + Y)]^ 1 , where w is so called division algorithm, that is 
w(&)w T (a) = a for any element a of the cone. The drawback of their extension was that the additional 
strong assumptions of invariance of t he distribution of the "quotient" by a group of automorphisms was 
imposed. To avoid this assumption Bobecka and Wesolowskil 2002 1 have developed another approach 
based on densities of X and Y. Assuming existence of strictly positive, twice differentiable densities on 
the cone of positive definite symmetric matrices they proved a characterization of Wishart distribution 
for U = (X + Y)~ X / 2 X{X + Y)~ x / 2 , where a 1 / 2 denotes the unique symmetric root of a positive definite 
matrix a. Among all division alg orithms the choice w (a) — a 1 / 2 is possibly the most natural one. 



Exploiting the same approach, lHassairi et al.l [20081 ] with the same technical assumptions on densities 



proved that the independence of X + Y and the quotient defined through the Cholesky decomposition, 
i.e. U = [W(X + Y)]~ 1 X[W T (X + Y)]' 1 , where W(X + Y) is an upper triangular matrix in the 
decomposition of X + Y, characterizes a wider family of distributions calle d Riesz- Wishart. This fa ct 
shows that the invariance property assumed in Olkin and Rubinl |l962 | and Casalis and Letad 19961 is 
not of the technical nature only. Analogous results for homogeneous cones were obtained in iBoutourial 



2009] and lBoutouria et al.1 [2011 1 



In this paper we attempt to lessen the smoothness assumptions in the IBobecka and Wesolows ki 2002] 
version of the Lukacs theorem. Our approach is directly through the Olkin-Baker equation and is very 
different from theirs. We succeeded in proving the characteriz ation theorem assuming continuity of 
densities only. It was possible thanks to an approach developed in lMolnarl [2006} . He proposed a method 
of solving functional equation for functions of matrix arguments consisting of two steps. First, an equation 
is solved for mutually orthogonal idempotents and afterwards a beautiful connection with Gleason's 
theorem is exploited. 
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This paper is organized as follows. We start with basic definitions and theorems regarding analysis on 
symmetric cones and a short introduction to Gleason's theorem in the next section. The statement and 
main result are given in Section Section [3] is devoted to solving two functional equations and is the 
technical core of the paper. It must be stated, however, that we limit our considerations to symmetric 
cones of rank strictly greater than 2 and exclude the octonion cone, due to the use of Gleason's theorem 

2. Preliminaries 

In the following section we give a short int roduction to the theory of sym metric cones and qu antum 
logic. For further details refer, respectively, to Faraut and Korarivil [1994 and Dvurecenskii 1993 1. 



A Euclidean Jordan algebra is a Euclidean space E (endowed with scalar product denoted (x,y)) 
equipped with a bilinear mapping (product) 

E x E 3 (x,y) (->■ xy G E 

and a neutral element e in E such that for all x, y, z in E 

(1) xy = yx 

(2) x(x 2 y)=x 2 (xy) 

(3) xe = x 

(4) (x,yz) = (xy,z) 

For an element x G E let L(x) : E — > E be linear map defined by 

L(x)y = xy. 

For x in E we define P(x) by formula 

P(x) = 2L 2 (x) -L(x 2 ) . 

The map P: E h> End(E) is called the quadratic representation of E. 

An element x is said to be invertible if there exists an element y in E such that L(x)y = e. Then y 
is called the inverse of x and is denoted by y = x _1 . It can be shown that x is invertible if and only if 
P(x) is invertible and in this case (P(x)) 1 = P (x _1 ). 

Euclidean Jordan algebra E is said to be simple if it is not a Cartesian product of two Euclidean Jordan 
algebras with positive dimensions. Up to linear isomorphism there are only five kinds of Euclidean simple 
Jordan algebras. Let K denote either the real numbers R, the complex ones C, quaternions HI or the 
octonions O, write S r (K) for the space of (r, r) Hermitian matrices valued in K, endowed with the 
Euclidean structure (x, y) = Trace (x • y) and with the Jordan product 

(1) L(x)y= i(x-y + yx) 

where x • y denotes the ordinary product of matrices and y is the conjugate of y. Then S r (M), r > 1, 
S r (C), r > 2, S r (H), r > 2, and the exceptional £3(0) gives the list of the first four kinds of Euclidean 
simple Jordan algebras. Note that in this case 

(2) P(y)x = yx-y. 

The fifth kind is the Euclidean space R™ +1 , n > 2, with Jordan product 

(3) (x , xx,. ..,x n ) (2/0,2/1, • ■ • ,Vn) = {^2 XiVu x ° yi + y ° Xl > • ■ • ,x ° yn + y° x ™j ■ 

To each Euclidean simple Jordan algebra one can attach the set of Jordan squares 

V = {x G E: there exists y in E such that x = y 2 } . 

The interior V is a symmetric cone. Moreover V is irreducible, i.e. it is not the Cartesian product of 
two convex cones. One can prove that an open convex cone is symmetric and irreducible if and only if 



THE LUKACS-OLKIN-RUBIN THEOREM ON SYMMETRIC CONES THROUGH GLEASON'S THEOREM 3 



it is the cone V of some Euclidean simple Jordan algebra. Each simple Jordan algebra corresponds to a 
symmetric cone, hence there exist up to linear isomorphism also only five kinds of symmetric cones. The 
cone corresponding to the Euclidean Jordan algebra M n+1 equipped with Jordan product © is called the 
Lorentz cone. 

We will now introduce a very useful decomposition in E - so called spectral decomposition. An element 
c G E is said to be a primitive idempotent if cc = c and if c is not a sum of two nonnull idempotents. 
A complete system of primitive orthogonal idempotents is a set {ci, . . . , c r } such that 

r 

Cj = e and acj = <5y C; for 1 < i < j < r. 

i=i 

The size r of such system is a constant called the rank of E. Any element x of a Euclidean simple Jordan 
algebra can be written asx = 5Zi=i ^ iCi f° r some complete system of primitive orthogonal idempotents 
{ci, . . . , c r }. The real numbers Ai, i — 1, . . . , r are the eigenvalues of x. One can then define trace and 
determinant of x by, respectively, tr (x) = yV., Ai and det(x) = ]X =1 Aj. An element x£E belongs to 
V if and only if all its eigenvalues are strictly positive. The rank r and dimV are connected through 
relation 

i. dr(r - 1) 
dimV = r + — y — -, 

where d is an integer called the Peirce constant. 

The Wishart distribution 7 Pi& in V is defined for any a 6 V and any p in the set 

{0, d/2, d,...,d(r- l)/2} U (d(r - l)/2, oo) 

by its Laplace transform 



exp(- (t, y})7 Pi a(efy) = (det (e + ta x )) 

v 



for any t + a 6 V. If p > dimV/r — 1 then 7 Pi& is absolutely continuous with respect to the Lebesgue 
measure and has density 

JpAdy) = ( d ^( a )) P ( det ( y)) p-dimV/ rexp( _ (a) y ))/v(y) rfy 5 

r v (p) 



where Ty is the Gamma function of symmetric cone V (see, for instance, Faraut and Koranvi |l994 |. p. 
122). 

One of the crucial elements in the proof of the main theorem is the use of Gleason's theorem which 
originated from quantum logic theory. Now we introduce necessary basics. 

Let L(H) be a set of all closed subspaces of a real, complex or left quaternionic Hilbert space H of 
finite dimension. We define charge as a mapping m: L(H) ->1U {— oo} U {oo} such that 

m(0) = 0, 




for a system of mutually orthogonal subspaces {at} teT from L{H) for any finite index set T. 

We denote by Pi(H) C L{H) the set of all one-dimensional subspaces of H. We say that charge m is 
Pi{H)-semibounded if 

inf {m{M): M e Pi(H)} > -oo. 
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Theorem 2.1 (Gleason's theorem) For every P\(H)-semibounded charge ni on L(H), 3 < dimiJ < oo, 
there exists a unique Hermitian operator t on H such that 

m(M) = Trace (t ■ p M ), M G L(H), 

where p M is an orthogonal projection onto M and Trace is the trace of a linear operator on H . 



For a f uller descript i on of above mentioned theory we refer to Dvurecenskii [l993j (real and complex 
case) and IVaradaraian 19851 ] (quaternionic case). 



3. Olkin-Baker Functional Equation 

In this section we solve the functional equation derived from the condition of independence of corre- 
sponding random variables. For details see Section [4] The following theorem is of independent interest 
in the functional equations theory and is the technical core of the paper. 

Theorem 3.1 (Olkin-Baker equation on symmetric cones) Let a, b, c and d be real continuous functions 
on a non octonion symmetric cone V of rank r ^ 2. Assume 

(4) a(x) + b(y) =c{x+y)+d(p({x+y)- 1 ^ as) , (x,y) eV 2 . 

Then there exist constants k\, k 2 G K, A G E, Cj G K, i G {1, 2, 3, 4} such that 

a(x) = (A, x) + ki logdet(a;) + C\, 
b(x) = (A, x) + k 2 logdet(a;) + C 2 , 
c(x) = (A, x) + (fci + k 2 ) logdet(a;) + C 3 , 
d(u) = k\ log det(u) + fc 2 log det(e — u) + C4 

{zeV:e-zeV}. 



(5) 



for all x G V and u G T> := 
The problem of solving 

(6) 



f(x)g(y) = p(x + y)q(x/y), (x, y) G (0, 00) 



for unkn own p o sitive functions f,g,p and q was first posed by lOlkinl 1975| . Its general s olution was 



given 



Baker 



1976 L and later, using a differ ent approach it was analyzed in lLajkol [l979j . Recently, 



n _ , n — — — -j . — j , — -i ^1 

Meszarosl [20101 ] and lLaiko and Meszarosl [20121 ] solved ((6|) assuming that it is satisfied almost everywhere 



on (0, oo) 2 for measurable functions, which are nonnegative on its domain or positive on some sets of 
positive Lebesgue measure, respectively. Finally, a new derivation of solution to when the equation 
holds almost everywh ere on (0, oo) 2 a nd no regularity assumptions on unknown positive functions are 
imposed, was given in iGer et ah 2011 1 . The equation ((4]) is a adaptation of © (after taking logarithm) 
to symmetric cone case. 

The proof is divided into two propositions and two lemmas. In Proposition 13.31 we show that each of 
the functions a, b and c is a sum of the function in the expected form, as in (|S|), and a homogeneous 
one, satisfying the same functional equation In Proposition 13.51 we show that these homogeneous 
functions are actually constant. 

We start with a simple lemma, which will be useful in the proof of the first proposition. 

Lemma 3.2 (Pexider functional equation on symmetric cones) Let a, b and c be measurable functions 
on a symmetric cone V satisfying 



(J) 



a(x) + b(y) = c(x+y), \/{x,y)eV 2 
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Then there exist constants a, /? € K and A G E smc/i i/ia£ V x G V 

a (a;) = (A, x) + a, 

(8) b{x) = (A, a;} + (3, 

c{x) = (A, a;) + a + p. 

Proof. First we will show that function c satisfies Jensen functional equation. For x, y G V we have 

a(x) + 6(y) = c(x + y) = a(y) + 6(x). 

Hence, a(x) — fe(x) = a(y) — 6(y) = const :— A\. Plugging a(x) = fe(x) + A\ into (J7JI for y = x G V we 
get 2fe(x) = c(2x) — A\ and in consequence 

c(2x) + c(2y) = 2c(x + y), 



i.e. fu nction c is Jensen on symmetric cone V. Following standard approach (see, for instance. iGer and Kominek 
201 lj ) we get that there exists a constant A 2 such that c(x) = /(x) + A 2 , where / is additive on V, i.e. 



/(x) + /(y) =/(x + y), x,yGV. 
We define an extension / of / to the whole E as follows 

/(x) for x G V 



/(x) = 



/(x + i x e) -/(t x e) forx^ V 



where i x = 1 — min^ Ai , A; being the i-th eigenvalue of x. Observe that for x ^ V element x + t x e belongs 
to V. Indeed, x + i x e has all its eigenvalues positive, actually greater than 1. Note that if x ^ V then at 
least one of its eigenvalues is non positive. Since t x > we have t x e G V. 

It can be easily verified that / is additive on E. Take for example x G V, y ^ V and suppose x + y ^ V. 
Then we have 

L = /(x) + /(y) = /(x) + f(y + t y e) - f(t y e) 

and 

R = /(x + y) = /(x + y + < x+y e) - /(< x+y e) 

Since 

/(x) + /(y + tye) + f(t x+y e) = /(x + y + t y e + t x+y e) = /(x + y + < x+y e) + f(t y e) 
we get that L — R. Hence, we know that there exists a constant A G E such that 

/» = (A,x) 

for all x G E. Since /(x) = /(x) on V the proof is complete. One gets the form of (jHJ) for A\ = a — (3 
and A 2 = a + f3. □ 

Proposition 3.3 Let a, b, c andd be real continuous functions on a symmetric cone V of rank r. Assume 

(9) a(x) + b(y) =c(x+y) + d(F^x+y)- 1 / 2 ^) as) , {x,y)eV 2 . 

Then there exist constants k\ , k 2 G R, A G E such that for all x G V and u G T> 

a(x) — (A, a;) + fci logdet(x) + e(x), 

b{x) = (A, a;) + fc 2 logdet(x) + /(a;), 

c(x) = (A, x) + (fci + k 2 ) logdet(a;) + g(x), 

d(u) = fci logdet(u) + k 2 logdet(e — u) + h(u), 
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where e, f , g and h are continuous functions satisfying 

(10) e(x) + f(y) =g(x+y) + h (p ((as + y)- 1/2 ) x 

for (x,y) G V 2 and 

(11) e(sx) = e(x), f{sx) = f(x), g(sx) = g(x) 
for any s G (0, oo) and x G V. 



Proof. The following proof adapts the argument given in lGer et al.l [201 1| . where the analogous result on 
(0, oo) was analyzed, to the symmetric cone setting. 
For any s > and (x, y) G V 2 we get 

(12) a(sx) + 6(sy) = c(s(x + y)) + d (p ((x + y)" 1 / 2 ) x) . 



Substracting now ([9]) from (|12|) for any s > we arrive at the additive Pexider equation on symmetric 
cone V 

o,(x) + b s {y) = c s (x + y), (x,y)eV 2 , 

where a s , 6 S and c s are functions defined, respectively, by a s (x) := o(sx) — a(x), 6 s (x) := 6(sx) — 6(x) 
and c s (x) := c(sx) — c(x). 

Due to continuity of a, 6 and c and Lemma 13.21 it follows that for any s > there exist constants 
A(s) £ E, a(s) G R and /3(s) G K such that for any x G V 

a s (x) = (A(s),x) + a(s), 
6 s (x) = (A( S ),x)+/3(s), 
c s (x) = (A(s),x)+a(s) + /3(s). 

By the definition of a s and the above observation it follows that for any (s, t) G (0, oo) 2 and z G V 

a st (z) = a t (sz) + a s (z). 

Hence, we have 

(13) (X(st), z) + a(st) = (A(t), sz) + a(t) + (A(s), z) + a(s). 

Since (fT3|) holds for any z G V we see that a(st) = a(s) + a(t) for all (s,t) G (0, oo) 2 . That is a(s) = 
ki log s for s G (0, oo), where k\ is a real constant. 
On the other hand 

(14) (A(si), z) = (A(s), z) + (A(i), «z) = (A(<), z> + (A(s), iz) 

where the latter equality is a consequence of insensitivity of the left hand side of (fT4"]l to the interchange 
of s and t. Putting s = 2 and denoting A = A(2) we obtain 

(A(i),z) = (A,z) (t-1) 

for t > and z G V. 

It then follows that Vs G (0, oo) Vz G V 

(15) a s (z) = a(sz) — a(z) = (A, z) (s — 1) + fci log s. 
Note that det(ax) = a r det(x). Define a new function e through equation 



a 



(x) = e(x) + (A,x) + ^logdet(x). 
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By (|15|) we obtain e(sx) = e(x) for s > and x 6 V. 

Analogous derivation shows that there exist functions / and g such that /(sx) = /(x), g(sx) = g(x) for 
s > 0, x G V and 

6(x) = (A,x) + ^logdet(x) + /(x) 
c(x) = (A,x> + ^ logdet(x) + ,g(x). 

for x G V. 

The definition of function h completes the proof 

h(x) = d(x) - ^logdet(x) - ^logdet(e - x). 

In order to show ([TP]) the Cauchy theorem on the determinant of the product det(P(x)y) = dct(x) 2 det(y) 
is used. □ 

The aim of the following proposition is to show that functions e, /, g and h obtained in the previous 
proposition are actually constant. As mentioned before we limit our considerations to symmetric cones 
of rank r > 2 except octonion cone, which allows to use Gleason's theorem. 

We start with a crucial lemma about the form of orthogonally additive function on the set of idempo- 
tents. Since every symmetric cone of rank r > 2 is isomorphic to r x r Hermitian positive semidefinite 
matrices over K, we will identify trace tr with the trace of a linear operator Trace . 

Lemma 3.4 Let E = S r (K) where K is either R, C or H with Jordan product (JlJ and V be the cone of 
E. 

a) Suppose that for a continuous function f : E — > R the equality 

f(p) + f(q)=f(p+q) 

holds for all mutually orthogonal idempotents p and q in E. If rank r is greater than 2 then there 
exists a unique Hermitian operator t£E such that 

f(p) = Trace (t-p), VpeE. 

b) If for any a?, y G V 

(16) Trace (t ■ log x) + 2 Trace (f • log y) = Trace (t ■ log (¥(y)x)) 

then there exists a constant such that t = $e. 

Proof, a) / is orthogonally additive M-valued function on the set of all idempotents. The set of all 
idempotents on E is compact in the norm topology, hence we obtain that / is bounded. Let H be a 
r-dimensional Hilbert space over K. For any closed subspace M C H define 

m(M) = f(p M ), 

where p M is an orthogonal projection onto M. Note that for any M G L(H) there exists a unique 
idempotent p M G E and any idempotent p G E is an orthogonal projection on a subspace of H. Hence, 
there is a unique correspondence between L(H) and the set of idempotents of E. 

It can be easily verified that m is a charge on L{H). From Gleason's theorem we conclude that there 
exists a Hermitian operator t such that 

m(M) = /(p M ) = Trace (t • p M ), VM C L(H). 

b) By ([2]) the ri ght hand side o f equation (fl6|) can be written as Trace (t • log (y • x • y)). Such equation 



was considered bv lMolnarl 2006] [eq. (2)], who proved the assertion for real and complex positive definite 



matrices only. His proof can be rewritten in our case, hence we get the result for elements of V also. □ 
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Proposition 3.5 Let e, f, g and h be continuous functions on non octonion symmetric cone V of rank 
r 2. Assume that (|10j) and hold true. 

Then functions e, f , g and h are real constants. 

Proof. The case r = 1 is trivial, since e(x) = e(sx) = e(l) for s = 1/x. 

Assume that r > 2. Pick any idempotent (not necessarily primitive idempotent) p on E and denote 
p 1 - = e p. Put x = ap + p 1 - and y = /3p + p x for (a, /3) S (0, oo) 2 . Since 

(x + y)- 1 / 2 = ( (a + /3)p + 2p^)- 1 / 2 = ^=p + ^ p-L 

we obtain 

e(ap + p X ) + /(/3p + p x ) - g((a + p)p + 2p x ) + h (^p + |p x 

This is a one-dimensional version of the main equation Q in functions of variables a and /3, so we already 
know its complete solution: 

e(ap + p x ) = A(p)a + Ki(p)log a + d(p) 
/(ap + p x ) = A(p)a + « 2 (p)log a + C 2 (p) 
5(ap + 2p- L ) = A(p)a + (/ii(p) + K 2 (p))log a + C 3 (p) 
M7P + |P X ) = «i(p) log 7 + «2(p) log(l - 7) + C 4 (p) 
C 1 (p) + C 2 (p)=C 3 (p) + C4(p) 



(17) 



for a > 0, 7 G (0, 1) and any idempotent peE. 

Observe that utilizing the homogeneity of function e and the above result we arrive at e(ap + p ) = 
e(^p x + p) and thus 

A(p)a + /ci(p) log a + Cx(p) = A(p x ) i + Kl (p x ) log ± + d(p x ) 

for all a > 0. From this we conclude that A(p) = 0, Ki(p x ) = — Ki(p) and Ci(p x ) — Ci(p). Analogous 
properties can be proved for k 2 and C 2 . Moreover, putting a = 1 in the first equation of (|17p we get 

e(e) = Ci(p), 

for all idempotents p G E. Therefore C,-(p) = C,*, i — 1, 2 for any idempotent p. 
Putting a — 2 in the third equation of dTTl) gives 

ff (2e) = («i( P ) + «a(p)) log 2 + C 3 (p) =: C 3 

for all idempotents p 6 E. This results in 

C 3 (p) = C 3 -Mp) + K 2 (p))log2 

and in consequence 

giap + p^) =g(2ap + 2p ± ) = («i(p) + «a(p)) log a + C 3 . 
We will derive now the formula for h(xp + yp x ) for (x,y) 6 (0, l) 2 . From equation (fT0|) we have 

e(ap + p x ) + /(p + /3p x ) = 5 ((a + l)p + (/3 + l)p x ) + ft ( -J-p + -^-p 1 

\a + 1 p + 1 

for a, (3 > 0. Taking x = anc [ y = — j-j- one gets 

/i(a;p + t/p- 1 ) = Ki(p) log - + K 2 (p) log - — - + C x + C 2 - C 3 

v i-y 

for x,ye (0,1). 
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Summarizing the above calculations, formulas for functions e, /, g and h are simplified to 

e(ap + p- 1 ) = «i(p) log a + Ci 
f(ap + p- 1 ) = K 2 (p) log a + C 2 
giap + p- 1 -) = (Ki(p) + K 2 (p))log a + C 3 
h(xp + yp^) = «i(p) log | + K 2 (p) log ±e§ + C 4 

^ Ci + C*2 = C3 + C4 

for a > 0, x, y G (0, 1) and idempotent p S E. Note that the above solution is also valid for cones of rank 
r = 2. 

Next step is to consider two mutually orthogonal idempotents p and q such that p + q / e. Such 
choice is possible only if r > 3, because there must exist at least three nonnull idempotents. 
Now we will rewrite the main equation (fTO]) in variables v = x + y and u = P ((x + y) -1 ^ 2 ) x. 

e(P(v 1 / 2 )u) + /(P(v 1 / 2 )(e - u)) = g(v) + h(u). 

Put v = tq + q 1 - and u = (3(tp + p- 1 ) for t > and /3, fit £ (0, 1). Then we have 

x = P(v 1 / 2 )u = 0(i(p + q) + (p + q)- 1 ), 

y = P(v 1 / 2 )(e - u) = (1 - 13) (tq + £f p + (p + q) 
After some easy but tedious computations one gets 



/(y) = /(«q+^fp + (p + q) X 

( 19 ) l-Bt 

= («2(q) + «i(q) + «i(p) - «i(p + q))log i + K 2 (p)log + C 2 . 



Due to the non symmetry in (jT9J) we have also 

/(y) = («2(p) + Ki(q) + «i(p) - Ki(p + q))log * — ^ + K 2 (q)log i + C 2 . 

fort > 0and/3,/3i G (0,1). Hence, «i(q) + Ki(p) = «i(p + q) holds for mutually orthogonal idempotents 
p and q on E. By continuity of / we get continuity of k\. By Lemma [3.41 a) there exists self-adjoint 
linear operator tx G E such that 

Kx(p) = Trace (ti ■ p) 

holds for any idempotent p. Similarly for k 2 (just substitute u by e — u in the above considerations) 

K 2 (p) = Trace (t 2 ■ p). 

Note that Trace (tj) = 0, because of the previously proven property fcj(p) = — ki(p ) for i = 1,2. 
The aim is to show that tx = t 2 = 0. Our first step to show this will be proving the following 



(20) 



5Z i=1 A « p *J = Z)i=i Trace (ti ■ p *) log Ai + Ci ' 

/ (E- =1 A * p *) = Trace (*a ■ ^) log A < + C 2, 

2( =1 AlP ^) = Z)^! TraCC ^ + ta ) ' P ^ l0g Ai + C3 ' 

^ (5Z i= i 7iP * ) = E i=1 ( Trace ^ tl ' p *) log71 + Trace (* 2 ' Pi ) log ^ ~~ 7i ^ + C4 
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for any system of orthogonal idempotents {Pj} i=1 such that J2i=i Pi = k < r, Xi > and 7.; G (0, 1) 
for i = 1, . . . , k. It can be easily shown using induction on k: 

• (basis) for fe = 3 it was previously shown for functions e and /. For g and h it is obvious since 

ff(x)=0(2x) = e(x) + /(x)-fc(£e), 
1 ' ft(x) = e(x) + /(e-x)- 5 (x). 

• (inductive step) Suppose now the formulas (|2T)|) are correct for k = n < r. We will show its 
correctness for n + 1 . 

Put v = tp n +p^anu = /3tp 1 + Y%=2 A Pi + £(Pi H ^ P n -i) X for t > Q, ft, /3t, fa, i = 1, . . . ,n - 1 

different and all in (0, 1). Then we have the following analogue of (|lgp 

n-l 

P^/^u = /3t( Pl + p„) + ^ f3 lPl + /3( Pl + ■■■ + pJ-S 

i=2 
n-l 

P(v 1 /2)( e - u) = (1 - /3i) Pl + " ft)p< + i(l - /3)p„ + (1 - /3)( Pl + • ■ • + p n ) x . 

i=2 

Now the result follows from similar calculation as in (|T9l) and (l2Tj) . 

Note that if x = J2i=i ^iPi then logx = YH=i Pi fog ^i? where logarithm a of a given element b G E 
is defined by exp(a) = b. 

We see now that formulas for functions e, /, g and h can be rewritten in a simpler form 

e(x) = Trace (ti • logx) + C\ 
/(x) = Trace (t 2 • logx) + C 2 
g(x) = Trace ((ti + t 2 ) • logx) + C 3 
h(u) = Trace (ti ■ logu) + Trace (t 2 • log(e - u)) + C 4 

C\ + C 2 =6*3+6*4 

for x G V and u G T>. Plug it into (fTU|) for (rx, y), r > 0, x, y G V. Since log (rx) = log (re) + logx = 
log (r) e + log x for any x G V, we get 

Trace (ti • log (x)) + Trace (t 2 • logy) = Trace ((ti + t 2 ) • log(rx + y)) 

+ Trace (ti • log (p ((rx + y)~ 1/2 ) *)) + Trace (t 2 • log (p ({rx + y) _1/2 ) y 
Passing to the limit as r — » and performing necessary computations we obtain for x, y G V 
Trace (ti • log (x)) - Trace (t x • logy) = Trace (ti • log (p (y _1/2 ) *)) 

and finally by putting y~ 2 in place of y we arrive at 

Trace (ti • logx) + 2 Trace (ti • logy) = Trace (t x • log (P (y) x)) . 



for x, y G V. From Lemma kOI b) we obtain that there exists real constant di such that ti = die. Since 
Trace ti = we get that i3i — 0. Similar argument holds for t 2 and thus the proof is complete. □ 



Note that Proposition 13.31 holds true with weaker assumptions, when equation is satis fied only almost 
ever ywhere for measur able, rather than continuous functions (see, eg. iGer et al.l 201 1| for (0, 00) case 
and iKolodzieiekl [2010j] for Lorentz cone case). The proof of Proposition 13.31 is valid for all symmetric 
cones including the Lorentz cone, while the proof of Proposition 13.51 is limited to cones of rank not equal 
to 2 and the octonion cone case is not covered. 

Proof of Theorem \3.1{ Follows immediately from Propositions 13.31 and 13.51 □ 
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4. LUKACS-OLKIN-RUBIN THEOREM WITH DENSITIES ON SYMMETRIC CONES 



The main result of the paper is the followin g theorem. We lessen the smo othness assumptions on 
respective densities from twice differentiability ( Bobecka and Wesolowski 2002] ) to continuity on sym- 
metric cones. Relaxing the regularity assumption in multivariate case to measur ability or even removing 
the assumption of density, w ithout invoking the invariance of the quotient as in lOlkin and Rubinl |l962j 
and Casalis and Letad [l996l |. remains a challenge. 

Theorem 4.1 (Lukacs-Olkin- Rubin theorem with densities on symmetric cones) Let X and Y be inde- 
pendent rv 's valued in a non octonion symmetric cone V of rank r > 2 with strictly positive and continuous 
densities. Set V = X + Y and U — P ((X + Y") -1 / 2 ) X . If U and V are independent then there exist 
a G V and p\,P2 > dim V/r — 1 such that X ~ "f Pl . a and Y ~ 7p 2 ,a- 

Proof. Let ip : VxV^PxVbea bijection defined by 



V(x, y) = (u, v) = (P ((x + y)- 1 ' 2 ) x, x + y) 



We have (U, V) = ip(X, Y). Inverse mapping rp 1 is defined by 



(x,y) = i/j 1 (u,v) = 

We are looking for the Jacobian of the map ip~ 

du\ /dx N 
rfv / \dy. 



,1/2 



,1/2 



that is, the determinant of the linear map 
/dx/rfu dx/dv\ ( du 



\dy/du dy/dvj ydv 



We have 



J = 



P( v l/ 2 ) 
-P(v!/ 2 ) 



dx/dv 
Idy — dx/dv 



A/ 2 \ 



dx/rfv 
Id v 



= Det(P(v 1 / 2 )), 



where Det denotes the determinant in the space of endomorphisms on V. From iFaraut and Koranvi 
[1994 Proposition III. 4.2] we get 



Det 



Now we can find the joint density of (U, V). Since (X, Y) and (U, V) have independent components, the 
following identity holds almost everywhere with respect to Lebesgue measure 

(det(x + y)) dimV/ 7x(x)My) = f v (p ((x + y)- 1 ' 2 ) x) / v (x + y), 

where fxjyju and fv denote densities of X, Y, U and V, respectively. Since the respective densities 
are assumed to be continuous the above equation holds for every x, y 6 V. Taking logarithms of both 
sides of the above equation we get 



(22) 
where 



a(x) + b(y) = c(x + y) + d (p ((x + y)- 1 ^ x) 



a(x) = log /x(x), 
b(x) = log /y(x), 

c(x) = log /v(x) - 
d(u) = log f v (u) 



dim V 



log dct(x), 



for x £ V and u G T>. 
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The thesis follows now directly from Theorem 13.11 Thus there exist constants fci, G K, A G E, 
Ci G R, i G {1, 2, 3, 4} such that 

/x(x) = e Q(x) = e Cl e <A ' x) (det(x)) fe \ 
JV(x) = e'W = e C2 e< A < x >(det(x)) fe2 

for all x G V. Since fx and /y are densities it follows that a = -A £ V, fc, = p, - dim V/r > — 1 and 

e C i = idetWZi . = 1 2 - ^ 

rV(Pi) ' ' 

Acknowledgement. The author thanks J. Wesolowski for helpful comments and discussions. 
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